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Multiple Polylogarithms: A Brief Survey 

Douglas Bowman and David M. Bradley 



Abstract. We survey various results and conjectures concerning multiple 
polylogarithms and the multiple zeta function. Among the results, we an- 
nounce our resolution of several conjectures on multiple zeta values. We also 
provide a new integral representation for the general multiple polylogarithm, 
and develop a g-analogue of the shuffle product. 



1. Introduction 

In recent years, nested harmonic sums have attracted increasing attention in 
both the mathematics and physics communities. The sums occur within the con- 
text of knot theory and quantum field theory yet their rich structure offers much 
to fascinate theoreticians in such diverse areas as algebra, number theory and com- 
binatorics. Multiple polylogarithms generalize the aforementioned nested sums, 
as well as the Riemann zeta function and the classical polylogarithm, while still 
retaining many interesting properties. Their study has led to many beautiful yet 
unproven conjectures, including evaluations at arbitrary depth discovered with the 
use of recently developed integer relations-finding algorithms and high precision 
numerical computations in the thousands of digits. 

Multiple polylogarithms 1461 148| are multiply nested sums of the form 

k 

(1-1) U su ..., Sk ( Zl ,...,z k ):= If/// :/ • 

ni>-">n.fc>0 j=l 

where s\, . . . , s k and z\, . . . , z% are complex numbers suitably restricted so that the 
sum converges. Instances of multiple polylogarithms have occurred in several 
disparate fields, such as combinatorics (analysis of quad-trees |41U59| and of lattice 
reduction algorithms [37j), knot theory [2H ggl ESI ED E3, perturbative 
quantum field theory [3 1101 120L I22| and mirror symmetry |55j . There is also quite 
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sophisticated work relating polylogarithms and their generalizations to arithmetic 
and algebraic geometry, and to algebraic iC-theory EH EH EfU H3 EHJ 023 
1ST] . 

Figuring prominently are the nested sums in which each Zj = 1. These 

latter are now commonly referred to as multiple zeta values |12l 1131 II 71 1681 182| 

and are denoted by 

k 

(1-2) C(*i,...,*fc):= E II" 

ni>--->nfc>0 j=l 

The study of such sums goes back to Euler 3Q] , who showed that 

m — 2 

(1.3) 2C(m, 1) =mC(m+l)- ^ C(™. - j)C(j + 1), 2 < m e Z. 

j=i 

It can be shown that the sum l|1.2|) is absolutely convergent in the region 

r 

{( Sl , . . . , s fe ) G C fc : E K ( s j) > r for 1 < r < ^1- 

i=i 

(The condition given in Proposition 1 of 83 is insufficient to guarantee absolute 
convergence.) 

Define the depth of the multiple polylogarithm to be the number k of 

nested summations. A good deal of work on multiple polylogarithms, and more 
specifically multiple zeta values, has been motivated by the problem of determining 
which sums can be expressed (say polynomially with rational coefficients) in terms 
of other sums of lesser depth. Settling this question in complete generality is cur- 
rently beyond the reach of number theory. For example, proving the irrationality 
of expressions such as £(5, 3)/£(5)C(3) appears to be impossible with current tech- 
niques. Nevertheless, considerable progress has been made with regard to proving 
specific classes of reductions, even at arbitrary depth. The first nontrivial success 
at arbitrary depth was the settling [121 113j of Zagier's conjecture |82| 

27T 4 " 

(1.4) C(3,1,3,1,--.,3,1) = 4-"C(4,4,...,4)= < n E Z, 

v v ' » v ' (4n + 2j! 

In n 

in which the 2n and n beneath the underbraces in l|1.4|) denote the depth of the 
respective multiple zeta values. Subsequent work (see eg. |16L 1171 Il8] ) has largely 
focused on developing a suitable framework for dealing with ultimately periodic 
argument strings in general, and additionally sums of multiple zeta values whose 
set of argument strings is fixed by the action of certain subgroups of the group of 
permutations. 

It is instructive to trace the development of the subject and see for oneself 
how ad hoc techniques and considerations have in many cases evolved into more 
systematic methods of ongoing interest. In this connection, one might begin by 
citing the partial fractions technique of Euler 39 and Nielsen |67| , subsequently 
employed by many others eg. [TTJJ [D>J EH ESI Ell E3E ES], and which Ohno re- 
cently parlayed in his exceedingly clever proof of the cyclic sum formula |57L I69| . 
Techniques based on elementary integration formulas and identities for special func- 
tions tailored to specific examples eg. 1101 1381 166] have evolved [111 I13j into 
quite general, sophisticated, and powerful analytic methods |16L II 8) . The naive 
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approach of deriving elementary series transformation identities and solving the 
resulting systems of linear equations jl 51 172) . used to prove reducibility results of 
depth three or less, has been largely superceded (eg. by methods based on contour 
integration 40 ) and supplanted by considerations of the shuffle and stufhe |13| mul- 
tiplications, and relatedly the harmonic algebra and the algebra of quasi-symmetric 
functions [^311541157] . 

Computational issues — both numeric and symbolic — have also come into play. 
Relations satisfied by multiple polylogarithms, and multiple zeta values in partic- 
ular, can be exploited by symbolic computer algebra systems to prove reductions 
of small weight |65j . (Here the weight of the multiple zeta value (|1.2I) is simply the 
sum of the arguments si + - • - + Sfc.) Interest in high-precision, rapid computation of 
multiple zeta values 35l 136) (see also |13l §7.2]) has been stimulated by the ability 
to numerically hunt for or rule out identities (to a high degree of probability) with 
the aid of recently developed integer relations-finding algorithms |42l 1501 I62| . 

In addition to the as yet unsolved problem of classifying all possible relation- 
ships between multiple zeta values at positive integer arguments, one can also con- 
sider (|1.2(l as a function of the complex numbers Si,...,Sfc and consider ques- 
tions regarding analytic continuation, trivial zeros, and values at the non-positive 
integers. The analytic continuation of (|1.2|1 in the case k = 2 was established 
by Atkinson pj] via the Poisson summation formula, and later by Apostol and 
Vu |3] , who used the Euler-Maclaurin summation formula. Subsequently, Arakawa 
and Kaneko [4] proved that if S2, . . . , s& are fixed positive integers, then (|1.2(l can 
be meromorphically continued as a function of si to the whole complex si-plane. 
The analytic continuation of 1|1.2|1 as a function defined on C k was established by 
Akiyama, Egami and Tanigawa pQ using the Euler-Maclaurin summation formula. 
An independent approach due to Zhao |83| uses properties of Gelfand and Shilov's 
generalized functions [44] . Zhao also attempts a discussion of trivial zeros for k < 3. 
To our knowledge, no-one has yet determined the trivial zeros of (|1.2(l for general 
k. 

The issue of values of l|1.2H at the non-positive integers is subtle, since for k > 1 
the result will in general depend on the order in which the respective limits are 
taken. Thus, for example, if n is a non-negative integer, s(k,j) and S(k,j) denote 
the Stirling numbers of the first and second kind, respectively, and Bj denotes the 
jth Bernoulli number, then [2] 



lim lim 

si—*—n s 2 ^0 



■ lim C(si)---)Sfc) 



(-1) 



n+ 1 




(-l) k +ij\S(n+l, 3 ) 
k + j 



whereas 




C ( Sl ,..., s ,) = (-i)^ n , -— — £ 

1 ; ' j=l 



s(k,j)B n+J 
n + j 




12 



5 



We are unaware of any systematic treatment in the case of arbitrary non-positive 
integer arguments. 
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1.1. Notation. Let cri,...,^ € { — 1,1}. We will have occasion to discuss 
the particular sums of the form 

k 

(i-5) u su ..., Sk (<T 1 x,<j 2 , ■ ■ ■ ,<Jk) = E xU1 n n i S3 a ] h > 

ni>--->ri fe >0 j=l 

in which < x < 1 is real and s±, . . . , are positive integers with x = S\ = o\ = 1 
excluded for convergence. Accepted practice dictates that 11.5fl may be abbreviated 
by Cr( s ii ■ ■ ■ i s k) with a bar placed over Sj if and only if Uj = — 1. When x = 1, 
these are called Euler sums. Thus a multiple zeta value is an Euler sum with 
no alternations. We adopt the convention that £ x () := 1 when no arguments are 
present (k = 0). We also drop the subscript x when x — 1 since Ci( s i> • • • i s k) 
agrees with ()1.2J) when each Sj is bar-free. For example, 

c(2,D = E i 4 L E 1 - 

n—l m— 1 

It will be convenient to abbreviate strings of repeated arguments by using exponen- 
tiation to denote concatenation powers. Then the first two members of (|1.4|) may 
be written C({3, 1}™) = 4" Tl C({4} n ). 

Finally, as customary the Gaussian hypergeometric function and the logarith- 
mic derivative of the Euler gamma function are denoted by 

F(a,b-,c-,x) = f:x-T\ { ; + j f + j l and ^ = 

respectively. We also abbreviate the set of the first k positive integers {1,2,..., k} 
byN k . 



2. Stuffles 

For the sake of brevity and simplicity, we shall restrict the discussion in this 
section to multiple zeta values. For a discussion of the more general polylogarithmic 
case, see |13) . 

As Hoffman |57j observed, one can view multiple zeta values as values of a 
homomorphism on a commutative Q-algebra in two ways; the Q-algebra multi- 
plications have been referred to elsewhere |13| as "shuffle" and "stuffle." It is 
conjectured that all relations between multiple zeta values are a consequence of the 
collision of the two multiplications, provided one admits the divergent sums (|1.2|) 
with s\ — 1 (suitably renormalizcd) into the model. However, there seems little 
hope of proving this conjecture in the near future, and at present a wide variety 
of analytic, algebraic, and combinatorial techniques are used to prove identities for 
multiple zeta values. 

Stuffle relations, or more simply stuffles — see i)3.1l and below for a discussion 
of shuffles — arise when one multiplies two nested series of the form (|1.2|l and expands 
the product distributively. Thus if u and v are (ordered) lists of positive integers, 
then 

c(u)c(«)= E cw, 
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where u * v is the multiset defined by the recursion 

(2.1) su*tv = s(u* tv) Ut(su*v) U (s + t)(s*t), Ks,(eZ. 

In (|2.1|l it is to be understood that if M is a multiset of lists and a is an integer, 
then aM denotes the multiset of lists obtained by placing a at the front of each list 
in M. For example (s, t) * u = {(s, t, u), (s, t + u), (s, u, t),(s + u, t), (u, s, t)} and 
correspondingly ((s, t)((u) — ((s, t, u) + ((s , t + u) + ((s , u, t) + ((s + u, t) + ((u, s, t). 

Let /(|i«|, \v\) denote the number of lists in u * v. The recursive decomposi- 
tion 112. 1 [I shows that the generating function 

oo oo 

F(x,y) :=^2J2f( m ^ xm y n 

m— n— 

satisfies the functional equation F(x, y) = 1 + xF(x, y) + yF(x, y) + xyF(x, y). It 
follows that F{x, y) = (1 — x — y — xy)~ l and hence that 

fc=0 v x v 7 fc=0 \ / \ / 

One can also give a direct, combinatorial proof of l|2.2|l by considering how the 
indices interlace in the product of two nested series of the form 11.2fl . 

There are interesting connections between stuffles, polyominoes, and codes 
which we briefly indicate. To begin, note that a stuffle counted by f(m,n) can 
be viewed as a pair ((f>,ip) of order-preserving injections 

(f> : N m -> N r , il>: N n -» N T 

where r is chosen so that max(m, n) < r < m + n and the union of the images 
of <j) and ip is all of N r . One can associate to such a pair a sequence of integers 
b±, . . . , b m by defining a\ — 4>{V) — 1 and a,j = 4>(j) — <fi(j — 1) — 1 for 2 < j < m 



and then letting 



bj 



—aj if (f>(j) is in the image of tp, 



dj otherwise 



for each j € N m . Since (j) is order-preserving, aj > for each j € N m and 
jyjLi = 4>{m) — m<n. Conversely, given a sequence of integers b%, . . . ,b m sat- 
isfying Y^jLi < n ; the pair ((p, ip) is uniquely determined. Let p = \{j : bj < 0}\. 
We have r — m + n — p, (f)(1) = \bi\ + 1 and </>(j) = 0(j — 1) + \bj\ + 1 for 2 < j < m. 
Put ■(/'(j) = (j>(j) if 6j < 0. The remaining values of tp are determined by the re- 
quirement that it be an order-preserving injective map of N n to N r . Thus, there is 
a one-to-one correspondence between the stuffles counted by f(n,m) and the sets 
of integer lattice points whose cardinalities satisfy 



(2.3) 



i(h,...,b m )eZ m -.jr\bj\<n\ = i(b u ...,b n )eZ n :J2\b 3 \<m\ 

I- j=l J L j=l > 



the identity l|2.3|l holding in view of the obvious symmetry /(m, n) = f(n, m). 

Define an n-dimensional polyomino formed by adding m coats to a single-celled 
polyomino, where a coat consists of just enough cells to cover each previously ex- 
posed (n— l)-dimensional face. There is clearly a bijection between such polyominos 
and the second set of lattice points (|2.3p . The relationship is explored in greater 
detail in [ISj • 
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3. Integral Representations 

3.1. The Drinfeld Integral. There is also a representation for multiple zeta 
values in terms of an "iterated" (Drinfeld) integral due to Kontsevich j82j . For real 
< x < 1 and positive integers Si, . . . , s& with x = s\ = 1 excluded for convergence, 
we have 

r k f S3 ~ Y rffW\ df U) 

<«) «-» «»>-/n(n3rh% 



r j ' / 1 — / 



where the integral is over the simplex 

x > 4 X) > • • • > i£> > • • • > t[ k) >■■■> *W > 0, 
and is abbreviated by 

(3.2) / TTa 5 ^ ''H, a = dt/t, b = dt/(l-t). 

Making the simultaneous change of variable t t— > 1 — t at each level of integration 
and then reversing the order of integration makes transparent the duality identity 
for multiple zeta values: 

(3.3) C(ai + 2, {l} r \ ...,s k + 2, {l} r >) = C(r k + 2, {1} S <% . . . ,n + 2, {l}* 1 ), 

first conjectured in 1521 and proved in |82) . 

A related integral representation enabled Ohno |68l to prove the following 
beautiful generalization of 1)3.3(1 . Let 

S(pi,...,p„;m) := ^2 C(pi + ci,...,p n + Cn), 

CiH hCji— m 

where the sum is over all non-negative integers ci,...,c n which sum to m. As 
in (|3.3J) define the dual argument lists 

p:= ( S i + 2,{ir,..., Si + 2,{ip) 

and 

p' := (r i + 2,{ip...,r 1+ 2,{in. 

Then |68| S(p; m) — S(p'; m). When m = 0, Ohno's result reduces to l|3.3|l . Another 
interesting specialization is obtained by taking p = (k + 1) and m = n — k — 1; 
one then deduces Granville's theorem 49 , originally conjectured independently by 
Courtney Moen |52| and Michael Schmidt [64] : 

t(si,...,s k )=C(n), 

siH hsjb=n 

where the sum is over all positive integers si, . . . , Sfe which sum to n and si > 1. 

The iterated integral representation is also responsible for a second multiplica- 
tion rule satisfied by multiple zeta values. Suppose that x, y E R and fj : [y, x] — > R 
are integrable functions for j = 1, 2, . . . , n. It is customary to make the abbreviation 



(3.4) / l[n / WfA'jUllj. a r .= fj{t s )dt s , 

J ~ L x>tt>t 2 >—>t n >y 3- L 
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with the convention that (|3.4|1 is equal to 1 if n = regardless of the values of x 
and y. There is an alternative definition of iterated integrals which explains their 
name. For j = 1,2, ... ,n again define the 1-forms ay by ctj := fj(tj) dtj. Then put 

.< f x I f x fi(t\) f' 1 a-i ■ ■ ■ a n dt\ if n > 

(3.5) / aia 2 ---a„ := < y ^ 

V ' Jy \ 1 ifn = 0. 

Expanding out this second definition, it is easy to see that it coincides with the 
definiton as an integral over a simplex. Both definitions occur frequently in the 
literature. 

Clearly the product of two iterated integrals of the form 1)3.411 consists of a sum 
of iterated integrals involving all possible interlacings of the variables. Therefore, 
if we denote the set of all (m + n)\/m\n\ permutations a of the indices N m+n 
satisfying cr _1 (j) < cr _1 (fc) for all 1 < j < k < m and m + 1 < j < k < m + n by 
Shuff(m, n), then we have the self-evident formula 

/ r x m \ / f x m+n \ r x m + n 

( 3 - 6 ) / n°y / n e / n^W' 

j=l 7 VJ « j=m+l 7 aeShuS(m,n) J y j=l 

and so define the shuffle product lu by 

(m \ , m-\-n \ m+n 

II «i):= E II " 
J = l ' ^j=m+l 7 o-eShuff(m,n) j = l 

Thus, the sum is over all non-commutative products (counting multiplicity) of 
length to + n in which the relative orders of the factors in the products a\a^ • ■ ■ a m 
and a m +ia„4-2 • • • a m+n are preserved. The term "shuffle" is used because such 
permutations arise in riffle shuffling a deck of to + n cards cut into one pile of to 
cards and a second pile of n cards. 

The study of shuffles and iterated integrals was pioneered by Chen 501 131] 
and subsequently formalized by Ree \74\ . As with the case of stuffles, one can 
view an element of Shuff (to, n) as a pair of order-preserving injections (4>, ip) where 
now <j) : N m — > N m+n and ip : N n — > N m+n have disjoint images. One can 
then define a vector (ai, . . . ,a m ) of non-negative integers by a\ = 4>(1) — 1 and 
a j = <fi(j) ~ <fi(j — 1) — 1 f° r 2 < j < to. Since is order-preserving, aj > 
for each j € N m and X^JLi a i = 4>{ m ) — m < n. Conversely, if we have such a 
vector of non-negative integers, then 4>(1) = a\ + 1 and 4>{j) = <j)(j — 1) + aj + 1 
for 2 < j < m defines an order-preserving injection : N m — > N m+n , and hence 
a shuffle. Thus, there is a one-to-one correspondence between Shuff (to, n) and the 
sets of non-negative integer lattice points whose cardinalities satisfy 



(ai , . . . , a m ) € Z> : aj < n 



(ai, . . . ,a„) G Z> : V^aj < to 



the latter identity holding in light of the fact that Shuff (to, n) is clearly symmetric 
in to and n. A deeper study of the algebra and combinatorics of shuffles leads to 
an alternative proof of (|1.4|l and generalizations thereof; see ^SJ 

3.2. A New Integral Representation. In light of the usefulness of the var- 
ious integral representations, it may be of interest to give here a new integral 
representation for The new representation appears to embody properties of 
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both the Drinfeld and partition integrals of [13] . and therefore may be useful in 
proving certain results for multiple polylogarithms which have thus far withstood 
attacks based on traditional methods. The derivation employs MacMahon's Omega 
operator, which discards terms with non-positive exponents from formal Laurent 
series in Ai, . . . , Thus, in view of Hl.ljl . if < x\, . . . , x k < 1, we may write 

Li S i,...,sfc • • • j x k) 

= E "/ f-''AV:>" • A ° :=1 

j = l nj >0 
k 

3=1 



«n n 



i>« 3 >--->«i >o 

1 s ; 

<i — 1 ^ i 1 / m ■ — 1 

J_ i>««> >..•>«<? >o J 

""^W I j=l V r=l M r 7 I mi>->m fe >0 j = l u Sj 

where A(s) denotes the set of all integration variables satisfying 

1 > Vlf 1 > 4 J) > • ■ • > > 

for j = 1, 2, . . . , k. Since < y 3 - < 1 for each j = 1, 2, . . . , k implies 



E \{y? 3 = E • • • E y- + -+^ 2/ - + - +nfc 

mi>--->m fe >0 j'=l ni = l n fc = l 

2/i 2/12/2 2/12/2 • • • Vk 



Tit. 



l-2/i 1 - 2/12/2 1 - 2/12/2 • ■ ' Vk ' 



it follows that 

(3.8) Li slv .. iS) , (xi,...,x k ) 
where t(x) := x/(l — x). 



J M?) J= i { V m =l y r=l 



4. Generating Functions 

In many cases, generating functions provide the best means of stating reduc- 
tions involving one or more parameters. A specific example of this which also 
illustrates how knowledge of the subject has progressed is given first. We then out- 
line a systematic approach for tackling multiple zeta values with periodic argument 
lists, followed by additional examples to illustrate the richness of the theory. 
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4.1. Two-Parameter Symmetry. In connection with Euler's result 1(1.3(1 . 
Markett [g4] derived 

C(*,l,l) = g«(« + l)C(» + 2)-i(.-l)C(2)((,)-i|;C( g - n -l)C(n + 3) 

n=0 

s — 4 n 

(4.1) + _J^( fl -„-2) J2 C(«-m + 2)C(m + 2), 3 < s G Z, 

n=0 m— 

via elementary but intricate series manipulations and partial fraction identities. An 
equivalent formula is proved in [10] using elementary facts about the dilogarithm, 
the polygamma function and the higher derivatives of the Euler beta function. For 
larger values of n, the representation of £(s, {1}™) in terms of values of the Riemann 
zeta function becomes increasingly complicated. Nevertheless, there is an elegant 
generating function formulation which we restate here. 

Theorem 4.1 ( jllp. The bivariate formal power series identity 

oo oo 

(4.2) ^^, m +V +1 C( m + 2,{lD 

771—0 n— 

= 1 - exp { jr \ (x k + y k -(x + y) k ) ((k)X 

^ k=2 ' 

holds. 

Corollary 4.2. Let n and s be non-negative integers with s > 2. Then 
C(s, {l} n ) lies in the polynomial ring Q[£(2), £(3), . . . , £(s + n)] . 

By comparing coefficients of x s ~ 1 y n+1 on both sides of 1(4.2(1 . one sees that in 
fact, £(s, {1}™) is a rational linear combination of products of Riemann zeta values 
such that the sum of the arguments in each product is equal to s + n. Moreover, 
Euler's result l|1.3|l is an immediate consequence of comparing coefficients of x s ~ 1 y 2 . 
Similarly, Markett 's formula l|4.1|l can be obtained most easily by comparing coeffi- 
cients of x s ~ 1 y 3 . Finally, as the right hand side of (|4.2|l is evidently symmetric in x 
and y, the left hand side must also be. Thus Theorem 14. II implies the special case 
C(m + 2, {1}") = C(n + 2, {l} m ) of the duality formula JjOJl. It would be interesting 
to find a generating function formulation of duality at full strength. 

4.2. Periodic Argument Lists. Results such as (|1.4H and (|4.2|l suggest that 
one might profit from a more systematic study of multiple zeta values whose ar- 
gument lists form an ultimately periodic sequence. This is indeed the case; such a 
study forms the basis of some of our current work in progress |18) . 

4.2.1. Period One. The case of all identical arguments is quite well understood. 
Nevertheless, there are a few items of interest worth recording here, in particular a 
connection to the problem of determining the number of unordered factorizations 
of an integer. 

For 3?(s) > 1, equation 11.2(1 implies 

°° °° / t s \ 

(4-3) £ <fcs C(M fe ) = n 1 + 77 • 
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If in (|4.3|) we take s to be an even integer, say s — 2n where n is a positive integer, 
then we may rewrite H4.j-SJ) in the form 

oo n— 1 

(4.4) J2(- 1 ) kt2kn m 2n } k ) = II sinc(Tr^), 

fc=0 j=0 

where p — e 7 ™/™ and sine a; = smx/x for i / 0; sincO := 1. The identity (|4.4|) 
is one of many possible generalizations of Euler's formula for £(2n), and moreover 
shows that (({2n} k ) is a rational multiple of Tr 2kn . 

Differentiating both sides of (|4.3[) and equating coefficients yields the recurrence 



fe 

(4.5) fcC(W fc )=£(-l) i+1 CC?*)C({«} fc ~ J )i 0<fceZ, U(s)>l, 
i=i 

which is really just a special case of Newton's formula 

fe 

ke k = *52(-iy +1 Pjek-j, 0<fceZ, 

3=1 

relating the elementary symmetric functions and power sum symmetric functions 

fe 



ji>--->jk>0 m=l j>0 



Substituting l/j s for each indeterminate Xj yields e k = C({ s } fc ) an d Pk = C(^ s )- 

From (|4.5|l it follows that if fc is a positive integer and 5ft(s) > 1, then £({ s } fc ) 
lies in the polynomial ring Q[£(s), £(2s), . . . , £(fcs)]. In fact, there is an explicit 
formula for C({ s }'°) m terms of a sum over partitions of k. 

Definition 4.3. Let r be a non-negative integer and let a — (ai, a^, ■ ■ ■ ) be a 
non-negative integer partition of r. Let rrij — ^{i : ai — j} be the number of parts 
of size j, and put c a = Oi>i m jK~ j) m> ■ Furthermore, abbreviate r = X)j>i a j by 
\a\ and Y[j>iPaj by Pa- 
in view of the generic relationship 63 

± e k t k = exp { - ± t*p£\ = 53(- t ) 

fe=0 r=l ' a 

for 3?(s) > 1 we therefore have 

f: t k a{ S } k ) = exp { - e ( - i) ^ (rs)r | = ^h)'"^- 1 n c(oj*)» 

fe=0 ^ r=l ' a ctj>0 

i.e. 

(4-6) C({*} fc ) = (~l) fe e c ^ n c(^ s )- 

a|=fc cy>0 
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We note the following connection with factorisatio numerorum |51| . (See 
also |29L 1701 EE].) Let a be as in Definition 14.31 Define the unrestricted divi- 
sor function associated with the partition a by 

d a (m) = 22 1. 

IliM d," J =m 

For example d% t i is the ordinary divisor function, and d 2 (m) = 1 if to is a perfect 
square and zero otherwise. 

Proposition 4.4. LetTk(m) denote the number of unordered factorizations of 
to into k distinct factors. Then 

r k (m) = (-l) fe ]T c- l d a {m). 

|a|=fe 

PROOF. Observe that for 5R(s) > 1, 

k oo 

C(M fe )= E l[nJ s = J2rk(m)m- s . 

n 1 >--->7i k >0 j — 1 m— 1 

Now compare coefficients of m~ s in l|4.6(l . □ 

Example 4.5. Since C({ s } 2 ) = \C{ S ) ~ K( 2s )' we S et T ^( m ) = ^di,i{m) - 
^d2(m). In particular T2(12) = 3. 

4.2.2. Period Two and Beyond. In contrast with the situation in which all 
arguments are identical, much remains to be explored in the case of argument strings 
of period two and higher. In |13| and |16| differential equations were found to be a 
useful technique for analyzing the generating functions for period 2. We summarize 
here some results from |13| and |16j to indicate the richness and complexity of the 
resulting formula? arising from the solution of the associated fourth order differential 
equation. 

Definition 4.6. For < x < 1 and z € C, let 

Y\(x,z) :— F{z, — z; 1; x), 

Y 2 (x,z) := (l-x)F(l + z,l-z;2;l-x), 

G(z) := \{ij}{l + iz) + i>{l-iz)-i>{l + z)-i){l-z)}. 

Theorem 4.7 ( 13 ). Let Y\ be as in Definition ^ .6] Then for < x < 1 and 

\z\<l, 

OO 

(4.7) E^ 1 )" 24 " 4 "^ 3 ' 1 }") =Yi{x,z)Y 1 {x,iz)- 

n=0 

Theorem 4.8 ([TS]). Let Y\, Y 2 and G be as in Definition \4.b\ Then for 
< x < 1 and \z\ < 1, 

oo 

(4.8) Y,(-l) n z An +H n ( x (3,{l,3} n )=G(z)Y 1 (x,z)Y 1 (x,iz) 

Yi(x, iz)Y 2 (x, z) Yi(x, z)Y 2 (x,iz) 
4Y 1 {l,z) 4Y!(l,zz) 
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Note that 

proves H1.4JL Similarly (|4.8(1 proves 

n 

C(3,{1,3}™) = 4""VC(4fc + 3)C({4}"- fc ) 



fc=0 



2n 4k ( 1 



^(4* + 2)! V 4 



C(4n-4A: + 3), 



which escaped the extensive numerical and symbolic searches carried out in the 
preparation of |HU12lll3| . Differentiation of l|4.8|l followed by a delicate analysis of 
the asymptotic behaviour of the requisite hypergeometric functions at their singular 
points proves |16 | the reduction 

C(2,{1,3}*) 

n , k ^ 

= 4-" ^(-l) fc C({4}"-' £ ) (4fc + l)C(4fe + 2) - 4 VJC(4j - l)C(4fc - 4j + 3) I 

fe=0 ^ 3=1 ' 

conjectured in I13j . 

The proof of Theorem 14 . 81 hinges on showing that both sides of (|4.7|l and (|4.8|) 
are annihilated by the same fourth order differential operator. In |13| . computer 
algebra was used to establish this for (|4.7|) . At the time, a conceptual proof was 
unavailable. Subsequently the present authors (see | 16p found a conceptual proof 
of the following more general result, which is perhaps best understood in the context 
of work going back to Orr |71j and Clausen 34 on differential equations satisfied 
by a product of hypergeometric series. The result is shown in |17j to be closely 
related to the combinatorial "shuffle" approach outlined in and may be stated 
as follows. 

Lemma 4.9. Let K be a differential field of characteristic not equal to 2 and 
let D be a derivation on K. For each k € K, define a derivation '■= kD. Let 
t be a constant, and suppose that for some f,g,u,V £ K the differential equations 
(DfD g +t)u — and (DfD g —t)v = hold. Then uv is annhilated by the differential 
operator {D 2 f D 2 g + At 2 ). 

In particular, taking f(x) = 1 — x, g(x) = x and t — z 2 , given that Y\ and 
Y-2 satisfy (DfD g + z 2 )y = 0, Lemma 14.91 shows that each of the three linearly 
independent functions Y\{x, z)Y\{x,iz), Yi(x, iz)Yi{x, z), and Y\(x, z)Y2{x,iz) are 
annihilated by the operator D 2 D 2 + 4z 4 . That L(x, z) and S(x, z) are annihilated 
by the same operator follows easily from the integral representation (|3.1fl . whence 
Theorem 14.81 is proved. 

Since D 2 D 2 +4z 4 is a fourth order differential operator, one might legitimately 
ask in what context the fourth linearly independent solution Y^ix, zyY^ix, iz) arises. 
It turns out that due to the double logarithmic singularity arising from the product 
of the underlying hypergeometric functions at x = 1, it is easier to ascribe a meaning 
to this solution in the case of alternating sums (|1.5(l . Recalling the generating 
function 



(4.9) A(*):=yVc({T} n ) = II 1 

n=0 j=l ^ 

from 1111. we have 



(-iy z \ r(i/2) 



r(i + z/2)r(i/2-z/2) 
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Theorem 4.10 ([IS]). Let < x < 1, and \t\ < oo. Put z = (1 + i)t/2, 
s = (1 + x)/2, and let U(s,z) = Y\{s,z) — zY2(s,z), where Y\ and I2 are as in 
Definition ^. b\ Then, 
00 

(4.10) J2 

Theorem 14. 101 is a bivariate generalization of the conjecture jl II equation (14)] 
in the case x = 1, and may be viewed as an analytic extension of the purely 
combinatorial identity (|5.4|) below. 

In recent work |18| . the authors have greatly extended the differential equation 
approach. The authors have obtained results on more general generating functions 
which include not only multiple zeta values, but polylogarithmic and hyperloga- 
rithmic |58| values in general. In fact, from the point of view of iterated integrals, 
arbitrary forms may occur in the iterated integrals studied. The differential equa- 
tions are still present. The authors have classified various bases for the solutions 
of the differential equations, given matrices for change of basis, and found the 
explicit representations of the monodromy matrices of the associated differential 
equations. These results actually stand out with greater distinction in a more gen- 
eral setting. Taking arbitrary 1-forms on a manifold M, an explicit homomorphism 
is obtained from ni(M, xq) into GL„(C). This gives rise to a transport between 
the manifold M and its principle bundle constructed from the representation into 
GL„(C). Finally these results can be cast yet more generally in the setting of dif- 
ferentiable spaces. Our homomorphism is similar to the celebrated homomorphism 
of K. T. Chen |3()l 1311 1321 133j in that it is built out of a generating function of 
iterated integrals. The essential difference is that Chen's homomorphism maps into 
a formal Lie group, while our homomorphism maps into GL„(C). Will our homo- 
morphism give different information than Chen's? We are currently investigating 
the geometric implications of our work in this area. 



t 2n U{iA} n ) + t 2n+i Uh{iA} n ) 



U (s, — z)U(s, iz) 
A(-z)A(iz) ' 



5. Shuffles and Cyclic Insertion 

As in 65 (cf. also |12l \74\ ) let A be a finite set and let A* denote the free 
monoid generated by A. We regard A as an alphabet, and the elements of A* as 
words formed by concatenating any finite number of letters (repetitions permit- 
ted) from the alphabet A. By linearly extending the concatenation product to 
the set Q (A) of rational linear combinations of elements of A* , we obtain a non- 
commutative polynomial ring with the elements of A being indeterminates and with 
multiplicative identity 1 denoting the empty word. 

The shuffle product (|3.7[) is alternatively defined first on words by the recursion 



(5.1) 



VlU £ A*, 1 LU W = W LU 1 = W, 

Va, b E A, Vu, v e A*, au lu bv = a(u lu bv) + b(au lu v), 



and then extended linearly to Q{A). One checks that the shuffle product so defined 
is associative and commutative, and thus Q(^4) equipped with the shuffle product 
becomes a commutative Q-algebra, denoted SIiq [A] . Radford 73 has shown that 
Shq[j4] is isomorphic to the polynomial algebra Q[X] obtained by adjoining the 
transcendence basis L of Lyndon words to the field Q of rational numbers. 
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The recursive definition l|5.1|l has its analytical motivation in the formula for 
integration by parts — equivalently, the product rule for differentiation. Thus, if we 
put a — fit) dt, b ~ g(t) dt and 

F(x) := [ (auuJbv)= ( f f{t) f udt)( f g(t) f vdt 

Jy \Jy Jy J \Jy Jy 

then writing F(x) — J* F'(s)ds and applying the product rule for differentiation 
yields 

F(x) = f X (f(s) f u] ( f g[t) f v dt] ds 



11 



g(s)[ I /(<) / udt) I vds 



y 



[a(u lu bv) + b(au lu v)] . 

Alternatively, by viewing F as a function of y, we see that the recursion (|5.1f) could 
equally well have been stated as 



(5.2) 



J Viu G A*, 1 lu = w lu 1 = ro, 

1 Va, b e A, Vit, v E A*, ua lu vb =(mlu vb)a + (ua lu v)b. 

Of course, both definitions are equivalent to l|3.7|l . 

The combinatorial proof 12 of Zagier's conjecture (|1.4|l hinged on expressing 
the sum of the words comprising the shuffle product of (ab) n with (ab) m as a linear 
combination of basis subsums. In |17| a more comprehensive study of the shuffle 
algebra SIiq [a, b] is undertaken, and as a consequence correspondingly deeper results 
for multiple zeta values are obtained. To highlight the most interesting of these 
results, we first recall the following 

Definition 5.1 (H2P- For integers m > n > let S m , n denote the set of 
words occurring in the shuffle product (ab) n lu (ab) m ~ n in which the subword a 2 
appears exactly n times, and let T m n be the sum of the m!/(2n)!(m — 2n)l distinct 
words in S m ,n- For all other integer pairs (m, n) it is convenient to define T mj „ := 0. 

One then has 

Theorem 5.2 (|17p. Letx andy be commuting indeterminates, and letm be a 
non-negative integer. In the commutative polynomial ring (Shoja, b])[x, y] we have 
the shuffle convolution formula 

m LW2J 

(5.3) xk y m ' k [( a6 ) fc m {ab) m ~ k ] = J2 ( 4x y) n ( x + y) m ~ 2n T ™.n- 

k=0 n=0 

A special case of Theorem 15.21 implies the intriguing shuffle factorization due 
to Broadhurst, and which in turn implies (|1.4|l : 

(5.4) A (jZ-) w A(^^j=M(z)G(Sh Q [a,b])[[z}}, i 2 = -1, 
where 

oo oo 

A(z) := ^2{z 2 ab) n (l + za) and M(z) := ^(z 4 a 2 6 2 )"(l + za + z 2 a 2 + z 3 a 2 b). 

n=Q n=Q 
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The experts will recognize (|4.9H as the analytic version of A(z) above, in which 
a = —dt/(l + t) and b — dt/(l — t). Similarly for M(z) and the left hand size 
of (|4.1U|) when x = 1. 

In addition, Theorem 15.21 plavs a key role in a remarkable combinatorial gen- 
eralization of H1.4fl which we proceed to describe. Let S mtn be as in Definition l5.ll 
Note that each word in S myn has a unique representation 

n 

(5.5) {ab) m ° W(a 2 b){ab) m ^b{ab) m2k , 

k=l 

in which too, mi, . . . , mm are non-negative integers with sum to — 2n. Conversely, 
every ordered (2n + l)-tuple (too, mi, . . . , mm) of non- negative integers with sum 
to — 2n gives rise to a unique word in S m ^ n via (|5.5|l . Thus, a bijective correspon- 
dence if is established between the set S m ^ n and the set Cm+\{m — 2n) of ordered 
non-negative integer compositions of to — 2n with 2n + 1 parts. In view of the 
relationship (j3.1|l expressing multiple zeta values as iterated integrals, it therefore 
makes sense to define 



Z(s) := / ip(s), s E C 2 „+i(to - 2n), 
Jo 



where as in (|3.2|) . we now identify the abstract letters a and b with the differential 
f-forms dt/t and dt/(l — t), respectively. Thus, if s = (too, mi, . . . , mm)-, then 



„i n 

Z[s) = \ (ab) mo \\{a 2 b){ab) m ^b{ab) 
Jo k=i 



= C({2} mo , 3, {2} m \ 1, {2} m \ 3, {2}"^, 1, . . . , 3, {2}" 12 - 1 , 1, {2}™-), 

in which the argument string consisting of rnj consecutive twos is inserted after the 
jth element of the string {3, 1}" for each j = 0, 1, 2, . . . , 2n. It turns out |17 | that 



27r 2m /m + 1 



(5.6) V Z{s) = , 

17 ^ W 2m + 2)!V2n+l 

seC2„+i(m-2n) 

for all non-negative integers m and n with to > 2n. The proof uses Theorem [5 
at essentially full strength combined with some tricky generatingfunctionology. 
Observe that equation l|1.4|l is the special case of l|5.6|l in which m = 2n, since 
Z({0} 2 «+ 1 )=C({3,l} n ). 

A more compelling formulation of l|5.6|l can be given as follows. Again, let 
s = (mo , mi , . . . , mm ) and put 

2n 

e(s) := Z(s) + ^ z (' n ^j> m j+i> ■ ■ ■ ,m 2n ,m , . . . ,TOj-_i). 

3=1 

In other words, sum over all cyclic permutations of the argument list s. Then |17| 
(5-7) Yl e(s) = Z(m) x \Cm+i(m -2n)\= I 

seC 2 „ + i(m-2n) ^ >' ^ 

is an equivalent formulation of l|5.6|l . Here, we have used 



Z(m) = C({2D= (2m+1) , , 0<meZ, 
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which follows from 1)4. 4[l . The cyclic insertion conjecture |13j can be restated as 
the assertion that C(s) = Z(m) for all s€ C2n+i(ni~2n) and integers in > 2n > 0. 
Thus, H5.7|) reduces the problem to that of establishing the invariance of C(s) on 
C2n+i(?ri — 2n). It is likely that this remaining step can be accomplished using only 
the shuffle property of multiple zeta values in conjunction with (|4.4|l . 



6. Dimension Conjectures 

Broadhurst 21] has conjectures concerning the size of various bases (graded by 
weight and depth) for expressing multiple zeta values in terms of either irreducible 
multiple zeta values, or irreducible Euler sums, and also for expressing Euler sums 
in terms of irreducible Euler sums. The adjunction of additional differential forms 
appears to simplify the problem at each stage. Thus, if Z?(n, k) denotes the number 
of multiple zeta values of weight n and depth k in a minimal Q-basis for reducing 
all multiple zeta values to a Q-linear combination of products of basis multiple zeta 
values, it is conjectured that 

n>3fc>l V ' y ' 

However, if we allow Euler sums into the basis, letting M(n, k) denote the minimal 
number of Euler sums of weight n and depth k needed to reduce all multiple zeta 
values to basis Euler sums, then 

n>3fc>l v Jy y ' 

One can also consider the problem of reducing Euler sums in terms of basis Euler 
sums. Let E(n, k) denote the minimal number of Euler sums of weight n and depth 
k required to reduce all Euler sums to basis Euler sums. It is conjectured that 

nri( i --v) w) =i- T ^. 

n>3 fe>l 

Adjoining forms associated with sixth roots of unity to the set of possible differential 
forms yields the multiple Clausen values |14| . and here it is conjectured that the 
number P(n, k) of irreducible multiple Clausen values of weight n and depth k is 
generated by 

2 

r n k\P(n,k) L i _ X V 



n>l fc>o 



7. (/-Shuffles 

Here we consider a q-analogue of the shuffle algebra discussed in ^ Let A be a 
set, not necessarily finite, and let rj : A —>■ A be bijective. Now form the free monoid 
generated by A and call it A* as before. Extend the action of rj to A* in the obvious 
way so that rj becomes an automorphism of A*. Again regard A as an alphabet, 
and the elements of A* as words formed by concatenating any finite number of 
letters (repetitions permitted) from the alphabet A. By linearly extending the 
concatenation product to the set Q{A) of rational linear combinations of elements 
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of A* , we obtain a non-commutative polynomial ring with the elements of A being 
indeterminates and with multiplicative identity 1 denoting the empty word. It is 
clear that r\ now extends to an automorphism of Q(^4). 

A g-shuffle algebra is defined to be the ordered pair (Q(A), lu q ) where lu q 
is a commutative and associative bilinear operator on Q(A) satisfying the identity 



(7.1) 



Vw e A*. 1 LU q W — W LU q 1 = W, 

Va, b S A, \/u,v(EA*, au lu q bv — a(u lu q bv) + b(rj(au) lu q v). 



We denote a q-shuffle algebra over A by Shqq [A] . It will be observed that a q-shufflc 
algebra is a commutative Q-algebra. 

Our definition implies that there may be more than one way of writing the 
q-shuffle product of two words. For example, letting a, b £ A, it is easy to see that 
in Shqq [A] one has o lu, b = ab + br/a = ba + arjb. As the length of the words 
being multiplied increases the number of different expressions also grows. 

The motivation for our definition of Shqq [A] is not difficult to see. As the 
shuffle algebra is motivated by the property l|3.6l) of iterated integrals, one wants a 
similar identity to hold for iterated Jackson q- integrals |43| . Recall the definition 
of a Jackson (/-integral. For x > 0, let / : [0, x] — » R be Riemann integrable. The 
Jackson q-integral of / on [0, a;] is defined by 

(7.2) f f(t) d q t := V f(xq n ) xq n (l - q). 

J ° n>0 

Because for any < q < 1 the sum on the right hand side of (j7.2|l is a Riemann sum 
for J Q f(t) dt, it follows that the Jackson g-integral tends to the ordinary Riemann 
integral in the limit as q approaches 1. 

One defines iterated Jackson g-integrals in exactly the same way that ordinary 
iterated integrals are defined by (|3.5|l . To this end, for j = 1,2, ...,n let fj : 
[0,x] — > R and u>j :— fj(tj) d q tj. Then put 

(7.3) / wiw 2 '-w„ := If/ fj{tj)d q tj, t := x 
Jo j=l^° 

Jo /l(*l) Jo' W 2 • • • d q tl if 71 > 

1 if n = 0. 

Here the fact that the 1-forms on the right hand side of l|7.3|l are q-difference 1- 
forms implies that the integral on the left hand side of (|7.3(l is a q-iterated integral 
and not an ordinary iterated integral. 

Iterating the definition of the q-iterated integral, one finds that 



(7.4) r Wl ...« k = ]T h{xq n *)h{xq n *+ n *--)---f k {xq" 



r,nk^ hni - 

x gfcnfc + (fe-l)nfc-lH h2n 2 +"i^ _ qj k x k 

but this is not a very convenient expression. To simplify (|7.4|l it helps to make the 
following change of indices: 

k 

li :— rij, 1 < i < k. 

j — k — i+l 
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Then (|7.4j) reduces to 
(7.5) / v 1 ---u>k = 



(l-q) k xq h f{xq h )xq h f{xq l *)---xq l *f(xq l «). 
0<h<h<---<l k 

This in turn motivates the definition of (/-difference 1-forms by the equation 

(7.6) u>i = Ui(t) = fi(t) d q t := /<(t)t(l - q). 

Notice that with this definition, the g-integral reduces to a summation operator on 
these 1-forms: 

pX 

uj = ^2oj(t Q q 3 ), 

j>0 

which agrees with the original definition by virtue of <o = x. 

Now fix < q < 1 and define the Rogers |75] g-difference operator 77 acting on 
continuous functions / : [0, 00] — > R by the equation (77/) (x) :— f(xq). Also define 
the g-derivative in the usual way: 

(D q f)( X ) := 

x(i-q) 

The relevant fact about the (/-derivative is the following well-known property: 

(7.7) D q f f(t)d g t = f(x). 

Jo 

We are ready to give the motivation for the q-shuffle product. The alphabet A 
now consists of g-difference 1-forms, and the automorphism 77 is Rogers' g-difference 
operator. The action of 77 is extended to forms by the equation rju; = to(tq) = 
f(tq)tq(l — q). This of course defines a new form u>' — r/uj by ui' = g(t)t{\ — q), 
where g(t) — qf(tq). The alphabet A needs to be infinite to account for all the 
forms rfu f or j g z. The action of 77 is now extended to Q(A) in the obvious way. 
It clearly forms an automorphism of this algebra. We wish to define the g-shuffic 
product so that for u,v e A* the following equation is true: 

(7.8) 

To accomplish this, one applies the g-analogue of the argument given in fJS] for 
deriving the recursive definition of the shuffle product. Take a, b £ A and u, v G A* . 
Put a = f(t) d q t, b = g(t) d q t and 



Jo 



F{x):= / {auv±s q bv)=[ / /(*) / ud q t)[ / g(t) / vd q t 
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Writing F(x) — f£ F'(s) d q s where F'(s) — (D q F)(s), and applying the g-product 
rule for ^-differentiation ((D q fg) — (D q f)g + (rjf)(D q g)) yields 

F(x) = J ( f{s) J Q U )(J Q 9{t) J vd ^) d ^ s 

+ [ 9(s)( [ q fit) [ udgt) f vd qS 
Jo \Jo Jo J Jo 

X 

[a[u uj q bv) + b(rj(au) lu v)] , 

where the first equality follows from (|7.7|l and the product rule for D q . Hence the 
inductive definition of the g-shuffle product results. Notice that commutativity and 
associativity of lu q follow immediately from (|7.8() . 

We will conclude with a few examples of the g-shuffle product illustrating how 
several equivalent sums can arise from this product. Taking wi lu,j LO2UJ3 using the 
inductive definition gives (among several possibilities): 

= +(^2(^1^3) + W2W 3 (l)Wl). 

Writing ojj = ajj (x) , these equations translate into the easily verifiable generic series 
identities: 

^WiOV 1 ) Yl U2{xq l -)LU Z {xq h ) 
0<h 0<l 2 <l 3 

2 LJ 1 {xq h )LU2(xq h )u J3 {xq l3 )+ ^ ^(^^(x^+VsCV 3 ) 
0<h<h<h 0<l 2 <h<h 

+ ^2(W 2 V 3 (V 3 Vl(V 1+2 ) 
0<l 2 <l 3 <h 

Y ^i(W> 2 (V 2 wV 3 )+ ^2(W 2 )^i(V i+1 v 3 (V 3+1 ) 

0<h<h<h 0<l 2 <h<h 

+ ^2{xq l2 )^{xq h )^i{xq ll+l ). 
0<l 2 <l 3 <h 

Taking the g-shuffle product as acting on the non-commutative polynomials in 
forms, it follows that all these different expressions for the g-shuffle product tend to 
the ordinary shuffle product in the limit as q approaches 1. Further results about 
g-shuffle algebras and their combinatorics will be given elsewhere. 
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